Probing spin correlations with phonons in the strongly frustrated magnet ZnCr 2 04 
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The spin-lattice coupling plays an important role in strongly frustrated magnets. In ZnC^CU, 
an excellent realization of the Heisenberg antiferromagnet on the "pyrochlore" network, a lattice 
distortion relieves the geometrical frustration through a spin-Peierls-like phase transition at T c — 
12.5 K. Conversely, spin correlations strongly influence the elastic properties of a frustrated magnet. 
By using infrared spectroscopy and published data on magnetic specific heat, we demonstrate that 
the frequency of an optical phonon triplet in ZnC^CU tracks the nearest-neighbor spin correlations 
above T c . The splitting of the phonon triplet below T c provides a way to measure of the spin-Peierls 
order parameter. 

PACS numbers: 75.50.Ee, 78.30.Hv, 75.30.Et, 75.10.Hk 



Geometrically frustrated magnets can resist magnetic 
ordering and remain in a strongly correlated paramag- 
netic state well below the Curie- Weiss temperature 0cw 
P, Q, . Particularly strong frustration is found in the 
Heisenberg antiferromagnet on the "pyrochlore lattice", 
wherein magnetic sites form a lattice of corner-sharing 
tetrahedra. The dynamics of the magnet is determined 
by the spin Hamiltonian 
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with the interaction restricted to nearest-neighbor bonds 
(ij). Theoretical investigations indicate that classical 
spins on this lattice may not order down to zero tem- 
perature 0,OJ- Cubic spinel ZnCr 2 04 offers a good real- 
ization of this model. Observation of magnetic order at 
low temperatures 0, Q has been explained in terms of a 
spin-driven Jahn- Teller effect 0, |j| that is reminiscent of 
the spin-Peierls instability in spin chains. The coupling 
between spin and lattice degrees of freedom thus plays a 
major role in relieving the geometrical frustration. 

The B sites of ZnCr204 are occupied by magnetic ions 
Cr 3+ with spin S = 3/2 [Fig. □»]. A crystal field of a 
nearly cubic symmetry splits the five 3d orbitals of Cr 3+ 
into a t2 g triplet and an e g doublet. By Hund's rules, the 
three electrons of Cr 3+ have aligned spins (S — 3/2) and 
occupy all of the three ti g states. The lack of orbital de- 
generacy precludes the ordinary Jahn- Teller effect com- 
mon to spinels and ensures a very small spin anisotropy 
0- The shape of the t 2g orbitals, pointing towards the 
neighboring Cr 3+ ions, makes direct exchange the pri- 
mary mechanism for magnetic interactions |lCt 1 1 11] . As 
a result, interactions beyond nearest neighbors are neg- 
ligibly small, while the nearest-neighbor interactions are 
highly sensitive to the motion of Cr 3+ ions creating a 
strong spin-phonon coupling. 

The magnetic susceptibility of ZnCr204 follows the 




(a) 



(b) 



FIG. 1: (a) One (out of five) of the symmetry coordinates of 
the Fiu mode which modulates Cr-Cr exchange, arrows show 
displacements; (b) spin ordering expected for a tetragonally 
distorted tetrahedron in ZnCr204, arrows indicate spins. 



Curie-Weiss law at high temperatures with &cw = 
390 K, which gives the nearest-neighbor exchange cou- 
pling J = 4.5 meV Q. As the temperature is lowered 
below Ocw, the magnet gradually enters a paramagnetic 
state with strong correlations between spins but no mag- 
netic order: the spins remain liquid but their motions 
are highly coordinated At T c — 12.5 K the mag- 

net undergoes a first-order transition into a phase with 
antiferromagnetic order and a structural distortion 0. 
According to the theory 0, 0, , the phase transition 
is driven by local distortions of the tetrahedra that have 
the E symmetry in the language of the tetrahedral point 
group Td- This is consistent with the observation of a 
tetragonal distortion below T c Q. 

In this work we demonstrate that lattice vibrations can 
provide quantitative information about spin correlations. 
We have measured the frequencies of IR-active phonons 
and found, below the magnetoelastic phase transition, a 
large splitting of a phonon mode involving magnetic ions. 
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FIG. 2: (Color online) Top curve: Rlt — reflectivity 
spectrum of ZnCr2C>4 single crystal in the low-temperature 
phase; bottom curve: 2(Rlt — Rht); inset: Rlt and Rht 
(line+symbols) in a narrow range. 



From the magnitude of the splitting we have inferred the 
absolute value of the spin-Peierls order parameter (S,; • 
Sj — Sfe • S;), where {ij) and (kl) are nearest-neighbor 
bonds 

Experimental details. Powders of ZnCr204 were pre- 
pared in air using the standard solid state reaction 
method. Small single crystals of ZnCr204 were grown 
using these powders using chemical transport method in 
quartz tubes sealed with CI2 gas as a transport agent. 
The crystals were of regular habit with the (111) face 
clearly visible. The working surface of our sample is tilted 
slightly (~3.5°) from the (111) plane as determined by 
X-rays. Reflectivity measurements were performed using 
a Fourier-transform spectrometer in the frequency range 
from 100 to 5000 cm -1 . The temperature dependence 
from 6 to 300 K was measured using liquid helium in a 
continuous flow cryostat (sample in vacuum) with optical 
access windows. 

The room temperature phonon spectrum of cubic 
spinels is well studied. Our measurements at room tem- 
perature agree with previously published data for sin- 
gle crystals of ZnCr 2 04 Q. As expected on symmetry 
grounds, the spectrum consists of four phonon triplets of 
the Fi u symmetry (point group Oh, space group FdZm). 
Figure [3 shows the reflectivity spectrum Rut for the low- 
temperature (LT) phase averaged over the temperature 
range from 9 to 11 K. It also shows a magnified differ- 
ence of the spectra 2(i?LT — -Rht), where Rht is the re- 
flectivity spectrum averaged between 13 and 14 K. The 
most prominent difference is observed at 370 cm -1 : be- 
low T c this phonon splits into two components. In ad- 
dition, there is a slight modification of the vibrational 
mode at 186 cm -1 , which cannot be fit with two slightly 
split oscillators. Several new modes with small oscil- 



lator strengths are seen above 300 cm~' in the low- 
temperature IR spectra. Most likely, these are phonons 
with nonzero wave vectors in the cubic phase that be- 
come visible in the LT phase due to Brillouin zone fold- 
ing |15l Il6j . This is consistent with a recent observation 
of an enlarged structural unit cell below T c [17j. Ra- 
man spectra also contain a set of new phonons in the LT 
phase [Hi ]. 

To extract the temperature dependencies of all phonon 
parameters, we fit the reflectivity R = |i-j-^I| 2 using a 
model dielectric function 

S 10 1 

e(u;,T) = e 00 + V- ¥~. (2) 

where e^T) = e(oo,T) is the dielectric constant well 
above all phonons, j enumerates the phonons, Sj(T) is 
an oscillator strength, uij (T) is a phonon frequency, and 
7j(T) is a damping rate. 

Figure |21 presents temperature dependences of the 
phonon frequency, oscillator strength, and damping rate 
for the 370 cm' 1 phonons. The symbols show the best- 
fit parameters for the model dielectric function J2J with 
four phonon modes above T c and five phonons below. 
Upon cooling from room temperature, the resonance fre- 
quency of the 370 cm' 1 phonon hardens first, as do the 
other three IR-active modes. In contrast, it softens signif- 
icantly below 100 K [Fig.Eta)]. Just below T c = 12.5 K, 
this phonon splits into two modes with a frequency differ- 
ence of 11 cm' 1 . The total oscillator strength is approx- 
imately conserved. There is no polarization dependence 
of the reflectivity spectrum at 14 K. Some polarization 
dependence appears on the split phonon below T c as can 
be seen from data points in Fig.[3]at 7 K. A clear observa- 
tion of the expected polarization effects in the tetragonal 
phase was not found most likely because of multiple do- 
mains in the sample. No hysteresis effect is observed. 

Theoretical model. The splitting of the 370 cm' 1 
phonon is consistent with the lowering of the crystalline 
symmetry from cubic to tetragonal. What physical mech- 
anism is responsible for the magnitude of splitting? An- 
harmonic effects unrelated to magnetism could cause the 
splitting, however there are two arguments against this 
interpretation. First, the magnitude of the anharmonic 
effect can be estimated as Alu/lo = -fAa/a, where 7 is 
an appropriate Griineisen parameter. To reconcile the 
observed splitting Aw/u = 0.030 with the magnitude 
of the tetragonal distortion Aa/a w 10~ 3 0, we must 
assume a rather large Griineisen constant 7 ~ 30, which 
must be explained. Second, the three other IR-active F\ u 
modes observed in the experiment exhibit much smaller 
(Aui < 0.2 cm -1 ) frequency related changes in the tetrag- 
onal phase. These objections cast doubts on anharmonic 
(and nonmagnetic) origins of the splitting. 

Most likely, the splitting of the triplet is caused by 
the same effect that triggers the spin-Peierls instability 
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- the spin-phonon coupling. It is well known that the 
elastic constants in magnetic materials are affected by the 
spins 0,H3- As a consequence, the phonon frequency is 
sensitive to correlations of spins of nearest-neighbor pairs 
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ujq + A(Sj • Sj), 



(3) 



where the constant A has a typical value of a few cm -1 
when optical phonon frequency is well above all magnons 
and phonon is not a zone folding mode; in addition to 
taking a thermal average, the spin correlations must also 
be averaged over the crystal with the weights appropriate 
for a given phonon. 

A strong argument in favor of the magnetoelastic 
mechanism is the nature of the 370 cm -1 phonon: it 
features by far the largest contribution of the symmetry 
coordinate modulating Cr-Cr bonds (Fig.Ek) among the 
four IR-active modes |2l| . The particular geometry of the 
occupied 3d orbitals in chromium makes the direct Cr- 
Cr exchange the primary source of magnetic interactions 
[ToL Hl|. In the case of direct exchange, the strength of 
magnetic interactions is particularly sensitive to the dis- 
tances between the magnetic ions (rather than the bond 
angles in the case of superexchange) . Therefore we ex- 
pect to find the largest splitting in those modes which 
produce the largest modulations of Cr-Cr distances, as 
indeed observed. 

To verify the linear relation Q and to determine the 
proportionality constant A, we have compared the tem- 
perature dependences of the spin correlations and phonon 
frequency in the high-temperature phase (Fig. EJ. The 
nearest- neighbor spin correlations (Si • Sj) were deter- 
mined from published data of specific heat Ne- 
glecting the magnetoelastic effects in the undistorted 
phase |23J we relate the spin correlations to the magnetic 
part of the specific heat per mole C m : 



(Si ■ Sj) = const + — l — [ C m (T) dT, 



(4) 



where 6Na is the number of Cr-Cr bonds per mole. The 
scaling relation works fairly well in the temperature 
range between 18 and 150 K. (The phonon softening 
above 100 K is likely due to thermal expansion.) This 
procedure yields the scaling constant A = 6.2 cm -1 . 

The relatively large value of A can be ascribed to the 
strong sensitivity of direct exchange to atomic displace- 
ments: the exchange constant J varies with the Cr-Cr 
distance approximately exponentially, 



J(r + Ar) « J(r)e" 



i Ar 



(5) 



on the scale of the Bohr radius a^. In contrast, the spin- 
phonon coupling in the case of superexchange comes from 
variations of bond angles, which have a less dramatic 
effect on magnetic energy. 
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FIG. 3: (Color online) Temperature dependence of the fit pa- 
rameters for the split phonon in ZnCr2 04; symbols — fit; 
squares, circles, and triangles at 7 K — fits of polarized spec- 
tra with P = 0°, 45°, 90°, respectively. 
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FIG. 4: Crosses: nearest-neighbor correlations (Si • Sj), up to 
an arbitrary additive constant, derived from the specific-heat 
data of Martinho et al. [2^|: straight line fits crosses below 
50 K. Squares: the phonon frequency. 



To obtain a quantitative estimate of A, consider a sim- 
plified lattice model in which magnetic ions of mass m are 
connected by springs of stiffness k (the rest of the atoms 
are discarded). The resulting "pyrochlore" lattice has 
only one optical mode that transforms as the irreducible 
triplet Fi u under the operations of the point group Oh- 
The spins move along a (110) direction [Fig. [Ja)] with 
the frequency luq — 2y/ko/m. The dependence of the 
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magnetic energy J(r)(S, ■ Sj) on the ion separation r 
provides a magnetic contribution to the spring constants 
k = k + Sk ex = k + J" (Si-Sj), where J" = d 2 J(r)/dr 2 . 
Hence A = 2J" /(mu ) « 2a 2 J/(mu) ). 

The value of a can be estimated theoretically from 
the exchange integral for 3d orbitals of hydrogcnlikc 
ions, which yields a — 2Z/3clb — 5.0 A for Cr 3+ 
ions (Z = 4) and A = 3.2 cm -1 . An alternative es- 
timate comes from the experimentally measured varia- 
tion of the exchange constant with the distance giving 
a = - J' I J = 8.9 A" 1 Ulllil and A = 10.0 cm" 1 . 
Our value A = 6.2 cm -1 is midway between these esti- 
mates. 

The observed splitting of the triplet phonon can now be 
translated directly into the spin-Peierls order parameter 
of the distorted phase. For a tetragonal distortion with 
lattice constants a = b ^ c (Fig. 

(Si • S 2 - S 2 • S 3 > = (w, - w»)/A = 1.8, (6) 

where bonds 12 are along the direction (110) and bonds 
23 are along (Oil). The bond variables (Si • Sj) are av- 
eraged both over the thermal ensemble and over the lo- 
cation in the crystal. 

Discussion. The inferred value of the spin-Peierls or- 
der parameter © is only a fraction of what can be at- 
tained in a state with a uniform tetragonal distortion and 
collinear spins (Fig. 5 of Ref.EI, in which bonds along 
the directions (110) and (110) are maximally frustrated: 
(Si • S2) — +S 2 , while the rest of the bonds are fully sat- 
isfied (Si -S 2 ) = — S 2 . In such a state, the z component of 
the phonon triplet [Fig. Ufa)] would probe the frustrated 
bonds with spring constants ko + J"S 2 , while the remain- 
ing two modes would involve bonds with k = ko — J" S 2 . 
We suspect that the primary reason for this significant 
reduction of the observed order parameter is the presence 
of a substantial nonuniform distortion in the crystal. Be- 
cause the spin-Peierls order parameter is linearly coupled 
to the crystal distortion , bond averages (Si • Sj ) will 
be nonuniform as well. Even though the local value of 
the order parameter © may be large, averaging over the 
crystal may reduce it substantially. 

Indirect evidence for a nonuniform distortion has been 
presented by Lee et al. Q. They point out that the ob- 
served uniform tetragonal distortion is too small to ac- 
count for the magnetic and elastic energy released in the 
form of latent heat at the phase transition. Direct evi- 
dence comes from X-ray measurements of Ueda et al. |J3 
and Raman scattering . If indeed the lion's share of 
the elastic and magnetic energy is hidden in the nonuni- 
form component of the lattice distortion, even larger 
splittings may be found in phonons with nonzero wave 
vectors. For a distortion with half-integer indices, such 
as (\\\) E3> the likely candidates are phonon doublets 
at the wave vectors 

In conclusion, we have studied the temperature evo- 
lution of optical phonons of symmetry Fi u in a strongly 



frustrated antiferromagnet ZnCr 2 04. The changes were 
most dramatic for the phonon with the largest am- 
plitudes of Cr-Cr vibrations. Its frequency showed a 
marked softening in the strongly correlated paramagnetic 
regime below 100 K. It essentially tracked the tempera- 
ture evolution of spin correlations inferred from magnetic 
specific heat allowing us to calibrate the frequency for the 
purposes of measuring spin correlations. Below T c the 
triply degenerate phonon splits into a singlet and a dou- 
blet providing direct information about the uniform part 
of the spin-Peierls order parameter 10. Its magnitude 
is only a fraction of what could be attained in a state 
with a uniform tetragonal distortion and collinear spins 
described in Ref. It is therefore likely that the lat- 
tice distortion and the spin-Peierls order parameter are 
nonuniform. The nonuniform component of the distor- 
tion may produce an even larger splitting of a phonon 
doublet with a wave vector {\\\). 
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